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Slow dynamics in glassy systems is often interpreted as due to thermally activated events between 
“metastable” states. This emphasizes the role of nonperturbative fluctuations, which is especially 
dramatic when these fluctuations destroy a putative phase transition predicted at the mean-field 
level. To gain insight into such hard problems, we consider the implementation of a generic back- 
and-forth process, between microscopic theory and observable behavior via effective theories, in a toy 
model that is simple enough to allow for a thorough investigation: the one-dimensional (p"^ theory at 
low temperature. We consider two ways of restricting the extent of the fluctuations, which both lead 
to a nonconvex effective potential (or free energy) : either through a finite-size system or by means of 
a running infrared cutoff within the nonperturbative Renormalization Group formalism. We discuss 
the physical insight one can get and the ways to treat strongly nonperturbative fluctuations in this 
context. 


I. INTRODUCTION 

Glass-forming liquids are systems whose salient phys¬ 
ical properties appear controlled by “nonperturbative” 
phenomena. At least in the deeply supercooled regime, 
when approaching the glass transition, the dynamics is 
best described as activated and cooperative.The pres¬ 
ence of such thermally activated processes is a proto¬ 
typical example where the role of fluctuations must be 
treated in a nonperturbative way. This is already well 
known in the case of standard homogeneous nucleation, 
for instance when a supersaturated vapor transforms into 
a stable liquid. The theoretical treatment of this prob¬ 
lem involves rare, localized events, which are described 
as nucleation droplets in a phenomenological approach^ 
and instantonic solutions of some free-energy functional 
at a field-theoretical level.® 

Activation is connected to “metastability”. In simple 
cases, the starting point of a theoretical description is a 
mean-held Landau free-energy functional, or a classical 
effective action in quantum held theory, which has sev¬ 
eral minima and is therefore nonconvex. This nonconvex¬ 
ity results from the absence of huctuations in the mean- 
held description and introducing huctuations in the the¬ 
ory leads to the exact free-energy or effective action with 
the needed convexity property. At the mean-held level, 
the deepest minimum is the stable thermodynamic state 
and the higher free-energy minima the metastable states 
(or false vacua in quantum held theory). The return to 
convexity has been theoretically described by explicitly 
accounting for excitations that are nonuniform in space, 
e.g. spin waves for systems with a continuous symmetry 
and droplets for Ising-like models, which encode the non¬ 
perturbative nature of the phenomenon. In the Ising-like 
case, nucleation (and growth) of a droplet then describes 
the escape of the system from a metastable state to reach 
the stable state and leads to activated dynamics.® 

In more involved situations encountered in systems 


undergoing a so-called huctuation-induced hrst-order 
transition,®’^ the mean-held theory does not provide a 
proper starting point as the relevant metastable states 
are themselves generated by huctuations. One must 
therefore hnd a way to include huctuations with wave¬ 
length up to some hnite scale in order to produce 
metastability and then study the localized excitations of 
such an effective theory.®’® 

The difficulty in the case of glass-forming systems is 
even stronger. The nature of metastability and of the 
metastable states is much more elusive®®^^^ and the ef¬ 
fective or coarse-grained landscape of minima and saddle- 
points is expected to be very complex, with a number of 
minima that is exponentially large in the system size. In 
this case, the “bottom-up” approach, deriving the be¬ 
havior of macroscopic observables starting from the mi¬ 
croscopic theory, i.e., interacting particles in the contin¬ 
uum, is just too difficult. Even computer simulations 
are of limited use for providing a full resolution because 
of the very fast growth of the equilibration time as one 
approaches the glass transition. 

A reasonable starting point would then be an effective 
theory that encodes the main physical ingredients while 
leaving out inessential ones. From a renormalization- 
group perspective, one would like to integrate out short 
length-scale fluctuations to obtain an effective theory 
that governs the long-range cooperative fluctuations. Nu¬ 
merical simulations can be particularly useful in this re¬ 
spect since they take into account, virtually exactly but 
for systems of limited size, all fluctuations. Provided one 
has some intuition about the nature of the effective the¬ 
ory and of the associated local order parameter, one can 
then try to extract the parameters of this theory from the 
exactly computed behavior of finite-size systems. This 
“top-down” approach would be instrumental in validat¬ 
ing and establishing the proper effective theory. 

Recently, there has been several numerical studies 
aimed at measuring the so-called Franz-Parisi potential 
U(g)^® in models of glass-forming liquids.This po- 
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tential plays the role of a Landau free-energy where the 
order parameter is taken as the similarity or overlap 
between configurations. It is at the root of recent ap¬ 
proaches that map the physics of supercooled liquids on 
effective theories of the random-field and random-bond 
Ising type.^®“^^ Studying for an ensemble of reference liq¬ 
uid configurations the average value of V(q) and its fluc¬ 
tuations in finite-size systems should provide a way to ac¬ 
cess some of the parameters of the effective magnetic-like 
theory, as discussed in Ref. 22. Once the parameters are 
known, predictions of the effective theory on long length 
scales can be further checked against other numerical or 
experimental observations. 

We think that this back-and-forth process, between 
microscopic theory and observable behavior via effective 
theories, is a key to solving the glass transition problem. 
It however requires a better understanding of the role 
played by nonperturbative fluctuations and the develop¬ 
ment of a theoretical approach able to capture them at 
all scales. The aim of this work is to study this problem 
on a toy model that is simple enough to be thoroughly 
investigated and use this as a benchmark for future work 
on glassy systems. 

The model that we focus on is the one-dimensional (p'^ 
scalar field theory defined by the Hamiltonian: 


HlW\ 





( 1 ) 


where L is the system size and the local potential V{(f) 
has a double-well form: 


= ( 2 ) 


We are interested in the low-temperature regime. There, 
the physics of the model is governed by strong nonper¬ 
turbative fluctuations: these are kinks or domain walls 
between the positive and negative phases. These local¬ 
ized defects have a finite cost and their density is always 
finite, albeit very small at low temperature (it follows a 
Boltzmann distribution). However, by redistributing the 
positions of these kinks the system can gain entropy. It 
is therefore their presence that destroys the phase tran¬ 
sition which is predicted at the mean-field level and re¬ 
mains in a perturbative treatment. 

Even though the present analysis is motivated by the 
simplicity of the model, which allows a detailed study, 
this is more than an academic problem. The one¬ 
dimensional Lp'^ field theory is actually central to several 
fields in statistical physics and condensed-matter physics 
where it appears mutatis mutandis in quite different 
problems: A Langevin dynamics in a double well,^^ quan¬ 
tum double wells^'^^^®, quantum-impurity problems^®’^*^ 
are all different incarnations of this very same model 
(with sometimes extra-difficulties and decorations). 

Let us now illustrate the main issues we are going to 
address in this work. 

The first issue concerns what we called the “top-down” 
approach. Mirroring the current situation in glasses, the 


problem we consider is one in which we want to infer 
the parameters of a theory that we conjecture to be of 
the Ising/i^^ type from the results of simulations and to 
further check that the effective microscopic theory we 
have in mind is the correct one. The input from simu¬ 
lations and other essentially exact computations that we 
consider as available knowledge is the probability of ob¬ 
serving a given average value </> of the field in a system 
of finite size L, or, more precisely, its logarithm, 

C/LW = -^ln/^L(</>) (3) 

where Pl{4') is the probability density to observe 
^ f dx(p(x) equal to ^ in a system of size L with peri¬ 
odic boundary conditions and j3 = l/{kBT). We will call 
Ul{ 4>) the finite-size effective potential since it takes into 
account all fluctuations exactly, up to the length-scale L. 
In the L —)■ 0 limit it coincides with the bare potential 
V{(f ), whereas in the thermodynamic limit it is equal to 
the free energy (exact effective potential) as a function 
of (f. This function ULif) is also known in quantum field 
theory as the “constraint effective potential”. 

Of course, in the case of the one-dimensional ip"^ field 
theory we know from the start that the proper effec¬ 
tive theory is just that given in Eq. (1). Nevertheless, 
the problem of inferring the bare parameters of the the¬ 
ory, c, r, and the energy of a kink from the behavior of 
finite-size systems is not straightforward. Moreover, un¬ 
derstanding in detail the evolution of Ul{4>), be., how 
the change in shape of the finite-size effective potential 
is related to the progressive integration of nonperturba¬ 
tive fluctuations is also very instructive. The knowledge 
gained in the case of this simple problem will likely be 
useful for tackling more difficult and still unsolved ones. 

The second issue is the development of a “bottom-up” 
approach that progressively takes into account fluctua¬ 
tions, including the nonperturbative ones, and allows one 
to eventually describe the macroscopic behavior. As ex¬ 
plained before, numerical simulations are not helpful in 
this respect since by construction they can be performed 
on finite-size systems only. This applies more specifically 
to glassy systems where the time scales needed to relax 
large systems close to the glass transition are unreachable 
even with the best available computers. Extrapolations 
to obtain the thermodynamic limit are then often dan¬ 
gerous and quite unrealistic. Needless to say that this is 
of course not true for the one-dimensional theory studied 
here. But as already stressed, the model is nonetheless 
used as a benchmark. 

The theoretical method of choice for progressively 
bridging the gap from microscopic to macroscopic physics 
is the renormalization group (RG).^^ Perturbative RG 
has been fully developed and understood since the 70’s 
and 80’s. The nonperturbative RG on the other hand 
has been the focus of an intense research only since the 
90’s: for reviews, see Refs. [9,33]. (At this point, we 
should acknowledge that there is always an ambiguity 
when using the adjectives “perturbative” and “nonper- 
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FIG. 1. Plot of Ul{4‘) as a function of (j> for different values of 
L, as obtained by MonteCarlo simulation: L = 4, 8,16, 32,128 
for model parameters r — 2, c = 4. 

turbative”. The former usually refers to an expansion 
in a few coupling constants and/or an expansion around 
the mean-field (gaussian) theory in powers of the differ¬ 
ence between the spatial dimension and the upper criti¬ 
cal dimension. The nonperturbative RG avoids such ex¬ 
pansions and is based on quite different approximation 
schemes that can potentially describe strong-coupling 
physics and, a key point for us here, the effect of non¬ 
perturbative fluctuations.^'*) 

The nonperturbative (NP) RG method that is cur¬ 
rently more used is the one introduced by Wetterich®’^®. 
It has been successfully applied to a variety of problems 
in high- and low-energy physics.In the field of statis¬ 
tical physics it has led to nontrivial solutions of long¬ 
standing problems in frustrated magnets,disordered 
systems such as the random-field Ising model,and 
out-of-equilibrium dynamical phenomena.'**’ The NPRG 
starts from an exact flow equation for the running ef¬ 
fective action, rfc[^!)], which is essentially the Legendre 
transform of the free-energy functional computed for an 
infinite system in which only fluctuations on length-scales 
less than 1/k have been integrated out. This running ef¬ 
fective action at (length) scale 1/k coincides in the ultra¬ 
violet or microscopic limit, fc —>■ A, with the bare Hamil¬ 
tonian and in the infrared or macroscopic limit, /c —>■ 0, 
with the exact effective action (Gibbs free energy) as a 
function of the field (/{x). 

The success of this NPRG relies on determining a 
simple yet rich enough truncation of the exact NPRG 
equation. There are cases in which this strategy has 
been instrumental in tackling problems in which non¬ 
perturbative fluctuations are present: the XY model 
in 2 dimensions,'** the random-field Ising model, 
or the return to convexity in the case of a first-order 
transition.However, the case we are considering in 
this work is somehow more difficult. There are strong 
nonlinear effects due to the interplay between the sharp 
changes in the field value on small length scales, which 
are associated with the kinks, and the long-wavelength 



FIG. 2. Plot of Uk{4>) = as a function of (j) for 

different values of k obtained by using the Local Potential 
Approximation (LPA). The model parameters are c = 1, 
/3 = 1, and r = 1/2. The flow equation is obtained with 
a regulator of the form*'* Rk{p) = {k^ — p^)0{k^ — p^). 
dkUk{4>) = —TT~^Uk {(f>)[Uk {4>) + Note that there is 

a fc-dependent but ^-independent contribution that is not in¬ 
cluded (hence the apparent difference of behavior with that 
in Fig. I). 

field variations on a scale of the order of the distance 
between the kinks. We find that this one-dimensional 
physics, which one of course knows how to solve by other 
techniques, is harder to access via the NPRG and remains 
an unsolved problem within this approach. We will com¬ 
ment in conclusion about the (more favorable) situation 
in higher dimensions. 

In both of the above situations, i.e., either in a system 
of finite size L or within the NPRG in a system in the 
thermodynamic limit but in the presence of an infrared 
cutoff on fluctuations of wavelength larger than 1/k, fluc¬ 
tuations are limited. As a result, the relevant potential, 
be it the finite-size one Ul{4>) or the rnnning effective 
one Uk{4>) = rfc[(/]/(/3L), need not be convex. Just like 
in the mean-field limit where no fluctuations are taken 
into account, which in the present case leads to a Landau 
potential equal to the bare metastability can thus 

be present. As the length scale over which fluctuations 
are allowed increases, i.e., with increasing L or decreas¬ 
ing k, metastability should become less pronounced and 
in the macroscopic limit, L —>■ oo or A: —>■ 0, both ULis/) 
and Uk{4>) should converge to the convex exact effective 
potential. 

The typical evolution with L of UL{(j>) is shown in Fig. I 
and that of the running effective potential Uk (</) is plot¬ 
ted in Fig. 2. (In the latter case, we have obtained the re¬ 
sult by using the so-called Local Potential Approximation 
(LPA)® of the exact NPRG equation.) The progressive 
disappearance of metastability is clearly observed in the 
two cases. The question we want to address in the former 
case is as follows: Say we are given some numerical data 
in the form of Fig. I; how can one extract information 
about the corresponding effective theory and its param- 
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eters? On the other hand, in the latter case we would 
like to develop an approximation to the NPRG that is 
able to reproduce the main features associated with the 
nonperturbative fluctuations in the present model, for in¬ 
stance the known fact that the curvature of the potential 
in 0 = 0 at low temperature is positive but very small as 
it behaves asymptotically as exp(—/3S'*), where S* is the 
energy cost of a kink. 

The rest of the paper is organized in two main sections, 
a first one where we address the top-down approach from 
finite-size studies and a second one where we discuss the 
bottom-up one through the NPRG. To avoid disrupting 
the flow of the presentation, some technical details are 
relegated to appendices. 

II. THE FINITE-SIZE EFFECTIVE POTENTIAL 

In this section we study the behavior of the finite-size 
effective potential and its evolution with the sys¬ 

tem size L. We first describe intuitively the shape of 
Ul{4') and explain the ideas on how to extract the rele¬ 
vant quantities such as the correlation length ^ and the 
surface tension 7 from the evolution of by focusing 

in particular on the behavior of two quantities: the cur¬ 
vature in (/) = 0 , = G£'( 0 ), and the height of the 

barrier between the potential in (/) = 0 and the minima in 
^ ~ ±1 (when present), = C/l(0) — C/l(± 1). We then 
present detailed analytic results for Ul{4‘) in the limit 
of zero temperature, which are obtained through the in- 
stanton technique, and use them as a benchmark to check 
the validity of our recipes for extracting the correlation 
length and the surface tension. Finally, we numerically 
determine the behavior of Ui,{(f) at finite temperature, 
i.e. at finite (but large) To this aim, we have com¬ 
bined Monte Carlo (MG) simulations and perturbation 
expansions based on a real-space RG and transfer ma¬ 
trix treatments. We then apply again the recipes for 
extracting the temperature dependence of ^ and 7 and 
compare the output with the direct numerical computa¬ 
tion of these quantities. 

A. The shape of Ul{4>) and its evolution with L 

At any given finite temperature, i.e., at any given finite 
correlation length if the system size goes to infinity, 
then the magnetization distribution goes to a Gaussian 
centered at (^ = 0 , due to the central limit theorem, and 
eventually converges to a Dirac delta function. (We use in 
this section the language of magnetic systems and call (j) 
the magnetization, or, more properly, the magnetization 
density.) Thus, in the thermodynamic limit, the finite- 
size effective potential displays a unique (parabolic) min¬ 
imum in (/) = 0. On the other hand, at any given finite 
system size, as the temperature goes to zero and the cor¬ 
relation length goes to infinity, the magnetization goes 
to either plus or minus one with probability one. For 



FIG. 3. Sketch of the evolution of the shape of the finite- 
size effective potential UL{<j>) for different system sizes L: a) 
L < 2cr; b) 2cr < L < G c) L < C; d) L ~ C; e) L > G f) 

L>e- 


^ ^ L the finite-size effective potential is given by two 
symmetric minima centered in ± 1 . Asa result, a nontriv¬ 
ial distribution of Tl( 0 ) and a nontrivial shape of U^icj)) 
arise between the two opposite limits considered above. 

In order to figure out intuitively the evolution of Pl {4>) 
and UL{(j)) with the system size, let us first focus on the 
typical configurations of the field (p{x) which, at least at 
low enough temperature, dominate the Gibbs measure. 
These are the configurations associated with the ground 
states of the system, corresponding to constant positive 
or negative magnetization profiles <f{x) = ± 1 , and the 
lowest excitations above them, involving domain walls 
{i.e., kinks and anti-kinks), which correspond to instan- 
tons that minimize the hamiltonian and connect posi¬ 
tively and negatively magnetized regions. At low enough 
temperature (large enough correlation length), the typi¬ 
cal configurations of the held are thus well described by 
regions with almost constant ±1 magnetization separated 
by narrow domain walls. The width tr of a domain wall is 
the typical size of an interface. The energy of a domain 
wall, S*, is, by definition, proportional to the microscopic 
surface tension 7 of the model, which is defined as the 
energy cost associated with the creation of an interface 
between two regions with opposite magnetization. It is 
easy to show (see below for more details) that the typi¬ 
cal distance between domain walls is of the order of the 
correlation length which is proportional to . 

If L is smaller than 2cr, no domain walls can be present 
in the system. (We consider periodic boundary condi¬ 
tions, so that the number of domain walls must be even.) 
Thus, Pl{4>) — andC/z,((/)) ~ V{(j)) (see Fig. 3a). 

Therefore, from Eq. (2), A ~ r/4 and k ~ —r. 

For L > 2a but still much smaller than the correla¬ 
tion length the probability of finding a domain wall 
is very small. The typical field configurations are then 
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approximately constant ±1 magnetization profiles plus 
some small thermal fluctuations, whose amplitude de¬ 
pends on V''{(j) = ±1) = 2r. On the other hand, con¬ 
figurations with zero magnetization correspond to field 
profiles with 2n domains walls, with n S N*, that are 
suitably placed between 0 and L. The thermodynamic 
weight of such configurations is proportional to 
and the probability of having (/) = 0 is obtained as 


Pl[4> = 0) oc 


{L-Aaf . 


g-2/3S* _|_ g-4/3S* 


2 L 

I {L - 6a)^ 

^ 192 


( 4 ) 


where the terms (L — 4cr)^/8, (L — 6tT)^/192, etc., cor¬ 
respond to the combinatorial factors accounting for the 
number of field configurations with 4, 6 , etc., domain 
walls between 0 and L that have zero magnetization 
(see the next section for more details). As long as 
2(7 <C T all configurations with more than a sin¬ 

gle kink/anti-kink pair are highly suppressed and their 
contribution can be neglected. Therefore, Tl((^ = 0) is 
dominated by field prohles with only two domains walls. 
Since all such profiles have the same combinatorial fac¬ 
tor (and thus the same probability), independently of the 
distance between the kink and the anti-kink, all interme¬ 
diate magnetization values, sufficiently away from ± 1 , 
occur with approximately the same probability. As a re¬ 
sult, for 2(7 <C A <C ^ the finite-size effective potential 
C/l(<(') is given by two deep narrow symmetric minima 
around ±1 (whose curvature is simply given by 2 r/A) 
that are separated by a central region where Ul {4>) is ap- 
proximatively constant. This is sketched in Fig. 3b. The 
barrier height = Ul{0) — Ul{H) is then given by 
2 S'*/L, and the curvature in ((> = 0 is = U’l{Q) ~ 0 . 

The qualitative shape of Ul (<(') does not show any sig- 
nihcant change until L < At this point, the terms of 
Eq. (4) corresponding to held conhgurations with more 
than two domain walls start to give a signihcant con¬ 
tribution to Pl{4>)- Since there are exponentially many 
more conhgurations of the domain walls corresponding to 
zero magnetization with respect to conhgurations yield¬ 
ing positive or negative magnetization, Al((^) starts to 
develop a secondary maximum around ()> = 0 as a result 
of this entropic effect. Correspondingly, develops 

a secondary minimum in zero, as sketched in Fig. 3c. In 
this regime the behavior of the barrier height and 
of the curvature kl are model-dependent and cannot be 
determined by simple intuitive argument: they must be 
computed in some explicit way, as we do in the following 
subsections. 

For L ~ ^ the barrier Al is expected to disappear as 
the value of the potential in (() = 0 crosses that in (^ ~ ±1 
(see Fig. 3d). As L further increases the minima in ±1 
become higher and eventually disappear. However, the 
potential may still remain nonconvex, as illustrated in 
Fig. 3e. Full convexity is recovered only for A —>■ oo. It is 
then easy to show that the hnite-size effective potential 
coincides with the Gibbs free-energy density (or exact 




FIG. 4. Schematic plot of the curvature kl = U'l{Q) (top) 
and of the barrier height Ai = Al( 0) — f7i(±l) times the 
system size A (bottom) as a function of A. The barrier height 
is shown on a log-log plot. The labels a)-f) correspond to the 
shapes of Ul{4’) of Elg- 3. Note that the behavior of AA_l 
at small A < 2 (t is not universal and depends on the bare 
parameters. 


effective potential) U ((/>) of the system: 

UL{,j,) = U{^) + o(^^y (5) 

where U{4>) is defined as the Legendre transform of the 
Helhmoltz free energy, 

u{ct>) = r^f{p,h) + H, ( 6 ) 

where h is the external magnetic held and {(p) = 

—df{l3,h)/d{l3h) = (j) (we have again used the magnetic 
language). As a consequence, for A ^ the hnite- 
size effective potential is a convex function of (j) and 
presents a unique minimum in (/) = 0: see Fig. 3f. In 
the thermodynamic limit the curvature kl approaches 
Koo = U"{Q) = where x is the magnetic susceptibil¬ 
ity dehned as x = d{'p)/d{ph)\h=o = L{{p^) - {(p)'^)\h=o- 
Based on the intuitive arguments discussed above, we 
can qualitatively determine the behavior of the quanti¬ 
ties of interest for us, kl and A^, as a function of A. 
They are schematically represented in Fig. 4. On very 
short length scales, A < 2cr, the curvature is negative, 
Kl — —r. Then, for 2(7 ^ A ^ kl is approximately 
zero. For L < ^, kl starts to grow and for A —>• oo it 
approaches 1/x as 1/A. In turn, the barrier height Al 
behaves roughly as 27 /A for 2cr <C A ^ ^ and rapidly 
vanishes for A > 
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We can therefore extract the important physical quan¬ 
tities by focusing on the behavior of kl and A/,. For in¬ 
stance, one possible recipe is to try to collapse the curves 
oi kl versus L obtained at different temperatures onto a 
master curve by rescaling the x and y axes by adjustable 
parameters. The parameters that provide the best col¬ 
lapse should then be x{T)~^ for kl {y axis) and ^(T) for 
L {x axis). Another possibility would be to plot L/S.^ as 
a function of L and, knowing the correlation length ^(T) 
from the previous operation, to look for a plateau or a re¬ 
gion of weak dependence on L for L < ^: at low enough 
temperature, the height of the plateau should then be 
twice the surface tension 7 (T). (Alternatively, one could 
do a log-log plot as in Fig. 4b.) In the next subsections 
we will implement and check these ideas in a quantitative 
way. 


B. Ul{4>) in the T —> 0 oo) limit 


As explained above, at very small temperature the 
Gibbs measure is dominated by the ground state of the 
system and the lowest excitations above it. The ground 
states of Eq. (1) correspond to constant field configu¬ 
rations <f{x) = ±1. The lowest excitations above the 
ground states correspond to nonuniform kink and anti¬ 
kink profiles that are obtained by minimizing the Hamil¬ 
tonian: 


SHl[(p{x)] 

6(p{x) 


V* 


= 0 


d‘^(p*{x) dV{ip{x)) 

dx"^ d(p{x) ’ 

(7) 


with the boundary conditions (p*{x —>■ —oo) = =Fl and 
(p*{x —>■ -boo) = ±1. This differential equation can be 
solved exactly for the Lp'^ theory in d = 1, yielding 


(p*{x) = ±tanh(x/CT), (8) 

with a = y^2cjr. The energy cost associated to such 
domain walls can be obtained by plugging Eq. (8) into 
Eq. (1), which gives S* = f HL[p>*{x)\dLX « y/8rc/9. 
Note that we can make more precise now the notion 
of “low-temperature regime”: It is obtained by let¬ 
ting either /3 or r be large, such that the Boltzmann 
factor associated with the presence of a domain wall, 
exp(—/3-\/8rc/9), is much smaller than 1. 

Small fluctuations of the field around the instantonic 
profile can be easily taken into account at a Gaussian 
level. After expanding the Hamiltonian around the in¬ 
stantonic solution up to second order in 6ip{x) = (p{x) — 
ip*{x), the thermodynamic weight of a single instanton is 
expressed as 


~ e 


-/3S* 


T>(p e 


-i / dxdy 


5<p{x)6<p{y) 


Scp(x) Sip{y) 


(9) 


In order to compute the functional integral above one 
thus need to diagonalize the operator corresponding to 


the kernel 


M{x,y) 


^ dx‘^ 


V"{^*{x)) 


5{x 


y), 


( 10 ) 


which yields 


Zi ~ e~^^* = e-/3'S*+5 ln(2’r/detM) _ 

In the present low-temperature limit, it can be shown 
that det M ^ j3r and, thus, S* ~ S* — (1/2)/?“^ ln(/?r) -b 

o(r'). 

At very low temperature the typical configurations of 
the field are therefore described by a dilute gas of domain 
walls separated by regions with constant (/? = ±1. The 
partition function of the system can thus be written as 
a sum over the number n of kink/anti-kink pairs (as dis¬ 
cussed above, the number of domain walls must be even 
to be compatible with the periodic boundary conditions) 
weighted by the energy cost times an appropriate 

combinatorial coefficient l 2 n accounting for all the pos¬ 
sible conhgurations of the positions of 2n domain walls 
between 0 and L: 


[L/2a] 

^L= ^2„(T)e-2"'5S*^ (12) 

n—0 

where [a;] denotes the integer part of x. Note that, since 
the instantons have a hnite width cr, we cannot place 
more than (Lj^a) kink/anti-kink pairs between 0 and 
L. The problem of determining the combinatorial coef¬ 
ficients l^n is equivalent to computing the entropy of a 
gas of 2n hard spheres of size ct on a ring of length L (see 
Appendix A). The resulting expression is 

72„(i) = ^ - 2na)2"-i. (13) 

n (2n — 1)! 

Two length scales thus naturally emerge from the cal¬ 
culation: CT, the typical size of an interface, and , 
which corresponds to the typical distance between two 
consecutive instantons and can be shown to be (twice) 
the correlation length ^ of the system. 

After introducing the rescaled variables C = Lje^^ 
and a = ct/L, the partition function finally reads 

[l/(2a)] 2n 

ZUC,t.) = 2 y: (14) 

The computation of the magnetization probability dis¬ 
tribution Pl{4’) in ths T —> 0 limit can be carried out in 
a similar way. Note that an analogous computation has 
already been done for the Ising model in d = 1^^ (see also 
below). 

For each given instantonic configuration with 2n al¬ 
ternate kinks and anti-kinks we define Xi, i = 1,..., 2n, 
as the lengths of the regions with constant <^9 = ±1. In 
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terms of these variables, the extensive magnetization M 
reads 

pL n 

M= ip{x)dx = ±'^{x 2 i-i - X 2 i) ■ (15) 

•^0 i=l 

Note that thanks to the translational invariance, one can 
choose without loss of generality to place the first domain 
wall at X = 0. The sign of M in front of the sum thus 
depends on whether the the first instanton is from ip = 
— 1 to v? = +1 or vice versa. Since each domain wall has 
a width a, we also have that 


2n 

Xi = L — 2na . (16) 

i=l 


In consequence, the extensive magnetization is bounded 
as \M\ < L — 2na. When enforcing the constraints given 
by Eqs. (15) and (16) one obtains 


Pl{M) = 


1 


SiM - L) + 6{M + L) 


l(L-\M\)/2a] 

+2 J2n(M,L)e-2"/5S* 

n—1 


(17) 


where is defined in Eq. (12). Again, the combina¬ 
torial factors J 2 n {M, L) can be computed exactly (see 
Appendix A). After introducing the rescaled variables 
a defined above and the magnetization density (jj = M/L 
and using the fact that S{L(j)) = {1/L)6{(j)) and Pl(L0) = 
(1/L)Pl((/)), we finally obtain: 


Pd^) = 


1 


ZLiC,a) 

[(l-|0|)/2a] 

+2 y] (c/2)2" 


- 1) + Si^ + 1) (18) 

[(1 — 2na)'^ — 


n!(n — 1)! 


where ZL{(^,a) is given in Eq. (14). It is easily checked 
that Pl((/') is properly normalized, Pl{(I}) dcj) = 1. 
One also finds that in the limit a —>■ 0, i.e., when the 
domain walls become infinitely sharp, and for S* = 2J, 
Eqs. (14) and (18) give back the exact results derived for 
the one-dimensional Ising model."*®. These calculations 
are explicitly done in Appendix A. 

The finite size effective potential and its evolution with 
the system size can be now explicitly determined in the 
T —>■ 0 limit from the relation in Eq. (3). t/L(^) be¬ 
haves as anticipated in the previous section: It presents 
two narrow minima in (j) = ±1, corresponding to the S- 
functions,^® and a secondary minimum in (^ = 0 due to 
the entropic term in Eq. (18). As L increases {i.e., a 
decreases) the minimum in (() = 0 becomes deeper and 
deeper, as the sum over n in Eq. (18) grows exponen¬ 
tially fast with C/ot. Eor L ~ ^ the value of the minimum 
in 0 = 0 crosses that of the two symmetric minima in 
(/) = ±1 (strictly speaking the value of the minima in 
(p = ±1 is defined only for a nonzero temperature;*® oth¬ 
erwise, one has to consider the weight of the delta peaks). 



m 


FIG. 5. Rescaled curvature as a function of L/^ for the 
ip'^ field theory in d = 1 in the low-temperature regime for 
P = 1, r varying from 2 to 6, and for c = 2r. 


Nevertheless, at any hnite L the potential remains non- 
convex due to the vestiges of the two minima in ±1. It is 
only in the thermodynamic limit that Ul^oo{4') = U{4’) 
recovers full convexity. 

From Eqs. (3), (14), and (18), one can compute all the 
desired characteristics of Pl (</>), such as the curvature in 
(j) = 0, kl, and the barrier height (when present), A^. 
Following the ideas presented in the previous section, we 
plot in Fig. 5 the curvature kl multiplied by the magnetic 
susceptibility x a function of the system size L divided 
by the correlation length for different temperatures.**" 
We have set (3 = 1, c = 2r and a range of r from 2 to 6 (as 
discussed above, the low-temperature limit here means 
that exp{—j3^/8rc/9) = exp(—4r/3) ^ 1). The curves 
show a perfect collapse, as expected. Via the instanton 
calculation we indeed have access to all physical quanti¬ 
ties of the present simple model. This is a consistency 
check of the recipe discussed above and an illustration of 
the range of temperatures where the asymptotic results 
apply. We show in Fig. 6 the evolution of the barrier 
height Ai multiplied by the system size L and divided 
by twice the domain-wall (free) energy S* as a function 
of for j3 = 1, c = 2r, and for r varying from 4 to 
10. This log-log plot is very similar to the sketch in Fig. 
4b. The value LAi^/{2S*) « 1 is observed for L ~ 2cr, 
which corresponds to very small values of £,/L, especially 
at low temperature. There is then a broad regime, up 
to p/L < 1, where one observes a small decay, by less 
than a factor of 10. Finally, for L > there is a fast 
(exponential) decay. 

These plots validate at a quantitative level the pro¬ 
posed ways of extracting the parameters of the theory 
from the behavior of the finite-size effective potential. 
We now turn to the same exercise but in the finite tem¬ 
perature regime where the analytical solution via the in- 
stantons is no longer a sufficient description. 
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FIG. 6. Log-log plot of the rescaled barrier, LAl/{2S*), as a 
function of for the field theory in d = 1 for /3 = 1, r 
varying from 4 to 10, and c = 2r. The curves stop on the low 
side for L = 2<j + a and the barrier is not defined within the 
instanton treatment for L < 2a. 


C. Ul{4>)/L for finite but iarge ^ 


In this section we apply and test the empirical recipes 
to extract x and 7 proposed above on a system at 
large but finite correlation length (corresponding to small 
but finite temperature). We obtain a numerical estimate 
of the finite-size effective potential Ul{4>) for a range of 
values of L through several methods. 

First, we have performed MC simulations. To this aim, 
we first discretize the continuum field theory of Eq. (1) 
by replacing the gradient by its discrete lattice version. 
The Hamiltonian thus becomes 


HLiWi}) = aJI ~ + 


(19) 


the instantaneous magnetization ip{t) = ( 1 /T) 
at regular time intervals corresponding to several times 
the equilibration time, say lOr. This allows us to make 
sure that the values of ipit) measured during the dynam¬ 
ics are statistically independent. In this way we con¬ 
struct a histogram of the magnetizations which gives an 
estimate of Pl{4’) and, from Eq. (3), we obtain UL{(j)). 
Results for /3 = 1, r = 2, c = 2r, and L varying from 4 
to 128 are shown in Fig. I. 

Note that in order to obtain an accurate enough es¬ 
timate of Pl{4‘) and of we need to sample rare 

events, which take place with an exponentially small 
probability in the system size. As a consequence, the 
number of measurements of the instantaneous magneti¬ 
zation must scale exponentially with L. Since the com¬ 
putational time of a single MC step scales linearly with 
the system size, this implies that the total computational 
time of our MC simulations scales as rLe^. Therefore, 
MC results are limited to not too large values of L, typ¬ 
ically L < 10^. 

In order to overcome this limitation and study larger 
system sizes, we have used a 1 /L perturbation expansion 
combined with an exact computation of the (Helmholtz) 
free-energy of the model through both a real-space RG 
approach and a transfer-matrix technique. 

Let us start with the definition of the magnetization 
probability distribution. 


where Tr{,^.j. = By using the integral 

sentation of the d-function, one easily obtains 

plOO 


( 20 ) 

repre- 


( 21 ) 


We set the lattice spacing o to 1 (note that for c = 2r 
the width of a domain wall is then a = 2 > a) and we 
consider periodic boundary conditions: ifiL+i = ‘Pi- 

The numerical simulations are performed with a 
Metropolis algorithm: At each time step we pick a site i 
at random, and attempt to change the value of (pi by a 
random quantity, 6ipi, extracted from a gaussian distribu¬ 
tion with zero mean and variance cr,^. We then compute 
the energy difference AHl = cSipi{S(pi + 2(pi — pi+i — 
Pi) + ^{Pi + ^Pi) ~ ^ {Pi) accept the move with the 

Metropolis probability p = minjl, Time is ad¬ 

vanced by 1/L. The typical width of the field shifts a^p is 
optimized recursively during the dynamics by enforcing 
that the acceptance rate of the moves (averaged over the 
last 100 MC steps) is approximately equal to 0.3. 

We start from a given initial condition (for instance 
Pi = -1-1 Vi) and let the system evolve and equilibrate. 
The equilibration time r, which of course depends on /3, 
r, and L, can be extracted from the exponential decay of 
dynamical correlation functions such as (1/L) '^i{pi{t + 
t')pi{t')) ~ {p^) e~*'P. In order to compute the mag¬ 
netization probability distribution, PL{(j)), we measure 


and 


-I plOO 

Ul{4>) = -T^ln / - /z,(/3,0), 

J — 200 

( 22 ) 

where fLiPip) is the Helmholtz free-energy density of a 
system of size L in the presence of an external uniform 
magnetic field /i//3: 

/l(/3,m) = . (23) 

For large enough L the integral in Eqs. (21) and (22) is 
dominated by the maximum in p = p*, which is given 
by48 


dPfL{P,p) 


dfi 


M* 


(24) 


Expanding the argument of the exponential around p* 
leads to 

(3fL{l3,p) + p(j) = /3fL{/3,p*) + p*(l)+ ^f?'\5pf 
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where /i"^ = 9 ”/l(/ 3,/i)/(9/x”|^* and Sfj, = fj, — One 
can thus treat all terms beyond the gaussian level in a 
perturbative way and obtain a systematic expansion of 
Pl{4>) and Ul{4>) in powers of 1/L. From a straightfor¬ 
ward calculation, one finds up to the order 1 /L^: 


UL{<l>)^fL{P,fin + J^-fL{P,0) 


(25) 


/ 27r 

1 


1 5[/F]2 ■ 


/3L2 

1 

OO 

' 24/3|/f^|3_ 


The above equation deserves some comments: 

( 1 ) As already mentioned, in the thermodynamic limit, 
Ul^oo{(P) converges to the Gibbs free-energy density 
U{(j>), which is defined as the Legendre transform of the 
Helmholtz free-energy density Jl^oo (/?, h) [see Eq. ( 6 )] 
and is therefore a convex function of the magnetization 


(2) Eq. (25) is actually an expansion in powers 
of ^/L. The successive derivatives of the Helmholtz 
free energy with respect to the external field /i yield 

(n) 

the n-points connected correlation functions, /?/* = 

which thus behave as 
As a result, the expansion of Eq. (25) does not 
converge for L/^ < 1 (even if L is large) and is expected 
to poorly behave compared to the numerical simulations 
in this regime. On the other hand, it should provide a 
good description of the finite-size effective potential for 

L/C > 1 . 

(3) In order to make some use of Eq. (25) we need to 
know the expression of the Helmholtz free energy of the 
model on a ring of L sites, at temperature /3, and in the 
presence of an external uniform magnetic field /r//3. 

The calculation of /i(/3,/r) can be done exactly by 
using a real-space RG approach, called the Migdal- 
Kadanoff (MK) scheme. It consists in integrating out 
iteratively half of the sites of the systems (say the odd 
sites) at each decimation step, and computing recursively 
the effective pair interaction potential, IF„(:/J, tp'), among 
the remaining sites. Consider for instance three consecu¬ 
tive sites, i, i + 1, and i -I- 2, at the p-th step of the renor¬ 
malization procedure. After integrating out the field on 
the site i + 1, one finds the following exact recursive equa¬ 
tion: 


IFp+l(pi,p*+ 2 ) = 


(26) 


^ J —oo 


with the initial condition: 

Woi(p,(p') = ^ {‘p-p'f 


(27) 


y{+) + y{+') - ^{+ + +') 


Eor a system of size L = 2^, after p— I decimation steps, 
there are only two sites left and the Helmholtz free-energy 
density can be obtained as a simple integration: 

-j ^+oo 

/l(/3,m) = In j_ . (28) 

This procedure allows one to obtain very accurate nu¬ 
merical values of fhifi, +) and of its derivatives, provided 
that the size of the system is an integer power of 2. In 
order to access other values of the system size L, we have 
complemented the RG calculation by a transfer-matrix 
(TM) approach. 

Indeed, the partition function of the system can be 
written as: 


p) — {ipi}'y(pi,v2'yvi,v3 ■ ■ ■ 


Vl,Vi 


= TrT^ = A[ -f A 2 + ... 


(29) 


where the transfer-matrix operator is such that = 

exp{—(3Wo{p,(p')) with Wq given by Eq. (27). One can 
then numerically diagonalize the operator by discretizing 
the values of the fields p and p' and compute its eigenval¬ 
ues, Ai, A 2 ,..., which leads to an approximate expression 
for the Helmholtz free-energy density. 


/l(/3,l) 



J_ „Lln(A2/Ai) 

I3l" 


+0 


(A 3 /A 1 



(30) 


Since the correlation length of the system is given by 


r' = -ln(A 2 /Ai) , (31) 


Eq. (30) provides a good approximation for /l(/ 3, p) only 
for L > 

The finite-size effective potential Ul{4>) is then ob¬ 
tained from Eq. (25). The numerical results for the cur¬ 
vature kl in (/ = 0 and for the barrier height at small 
but finite temperature (or rather, correlation length) are 
displayed in Eigs. 7 and 9. In Eig. 7a, we plot kl versus 
L for several temperatures (actually, values of r as we fix 
/3 = 1) and in Eig. 7b we show the best data collapse on 
a mastercurve after rescaling both the curvature and the 
system size by temperature-dependent adjustable param¬ 
eters. (Note that the curvature is obtained from the 
1 /L expansion only as the numerical accuracy of our MG 
data is not high enough to allow a good determination of 
the curvature.) In Eig. 8 we plot the best-fit parameter 
^fit versus r and compare it to a direct determination 
of the correlation length through MC simulations and 
the instanton technique: we find a very good agreement 
between the two sets of data. The same agreement is 
obtained for Xfa which is found proportional to the cor¬ 
relation length, ^fit or as expected in one dimension. 

In Fig. 9a, we display a log-log plot of LA^ versus 
L/^fit where is obtained from the previous data col¬ 
lapse in Fig. 7. (As could be anticipated, the 1/L expan¬ 
sion fails completely for L/^ < 1 and is not shown here.) 
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FIG. 7. (a): Log-log plot of the curvature kl = U'l{Q) as 
a function of L for /3 = 1, r = 4, 6, 8, 10, 12, and c = 
2r. (b): Same data with a rescaling of the x and y axes, as 
L/^fit and respectively, to provide the best collapse 

to a mastercurve. 


FIG. 9. (a): Log-log plot of the barrier height times the 

system size, LA.l, versus L/^fu for /3 = 1, c = 2r, r = 4, 
6, 8, 9, 10. Filled circles correspond to MG data. Empty 
squares are obtained using the 1/L expansion with the real- 
space RG approach and filled triangles correspond to the 1/L 
expansion with the transfer-matrix approach, (b): Same data 
with an adjustment of the y axis, LAl/xiu, to provide the 
best collapse for L/^fu < 1. 



FIG. 8. Plot of the scaling parameter ^fu versus r and com¬ 
parison with the direct computation of the correlation length 
through MG and instanton techniques. Recall that have set 
/? = 1, so that the temperature dependence is controlled by r. 
The error bars associated with the htting procedure are very 
small in this case and of the order of the symbol size. 



FIG. 10. Plot of the scaling parameter 7 /it, multiplied by P/r, 
versus r and comparison with the direct instanton computa¬ 
tion of the surface tension. The error bars are associated to 
the uncertainty of the collapse procedure. Recall that we have 
set /? = 1 and c = 2r, so that the temperature dependence 
is controlled by r, with 7 oc r. Note also that the instanton 
calculation breaks down at small r. 













11 


Fig. 9b then shows the same data with LAl divided 
by a temperature-dependent adjustable parameter 
that ensures the best collapse of all curves for L/^ < 1 
(this parameter is determined up to a multiplicative con¬ 
stant). When plotted as a function of r, we find that 
this best-fit parameter "ffa matches very well the r de¬ 
pendence of the direct estimate of the surface tension of 
the model through the instanton technique: see Fig. 10. 
Here, /3 = 1 and c = 2r, so that (33* = 2r/3 — (1/2) Inr 
for large enough r. We have arbitrarily adjusted the un¬ 
known constant in 'jfn so that the latter is roughly equal 
to 2(33*: the plot in Fig. 9b is shown with this choice 
of constant (which merely shifts all curves by a constant 
amount on the log scale). 

These plots therefore confirm that the ideas and recipes 
we have proposed to extract the correlation length the 
susceptibility x and the surface tension 7 (or alterna¬ 
tively the amplitude of the gradient term c) from finite- 
size numerical data for the effective potential work nicely. 
Without a priori knowledge one can empirically deter¬ 
mine the relevant parameters of the underlying (effective) 
theory from observations on finite-size systems. 


III. NONPERTURBATIVE RG 

In this section we take a different approach than the 
above one. It is more of a “bottom-up” approach where 
we start from a known microscopic (or to the least effec¬ 
tive) theory and try to include the fluctuations, in partic¬ 
ular possibly strongly nonperturbative ones, to describe 
the observed macroscopic behavior. As already stated, 
the method of choice to achieve this is the RG, more 
precisely the nonperturbative RG (NPRG). 

The NPRG emerged from Wilson’s work^^ in the early 
70’s and has since been formulated is several alternative 
approaches.All of them have the common denom¬ 
inator of treating and summing up fluctuations in a con¬ 
tinuous way. In this work we focus on the formalism that 
has been originally developed by Wetterich and cowork¬ 
ers since the 90’s.®’^® In a nutshell, the idea is to start 
with a given held-theoretical model described by a micro¬ 
scopic (bare) action and to add to it an infrared 

(IR) regulator in the form of a mass term: 

S'):/?] ^ 3k[:p] = S'):/?] + ^ (p{q)Rkiq)(p{-q) (32) 

where f = f dq/{2Tr)^ with D the space dimension and 
fc is a running (momentum) scale. The IR cutoff function 
Rk{q) goes to zero when fc —>■ 0 and provides a mass to 
the small-momentum modes, Rk{q) ^ bkk^ when g —>■ 0 , 
but is otherwise arbitrary. One also has that k < A where 
^ 1 is the shortest wavelength on which the field (p 
can fluctuate; A is the ultraviolet (UV) cutoff scale where 
the continuum theory meets the microscopic details. 

From the regularized action 3k [p] one can define a cut¬ 
off dependent generating functional of the connected cor¬ 


relation function (the analog of a Helmholtz free-energy 
functional at the scale fc). 


ITfc[J] = Iny exp -(33k[p\-\- J J{x)(p{x) 
where f^ = f d^x, and its Legendre transform. 


, (33) 


FfcM -b Wk[J] = J J{x)(j){x) - ^ j (l>iq)Rk{q)((){-d), 

(34) 

where for convenience one subtracts the contribution 
from the regulator in the definition of F^ [(/]. The latter is 
called the effective average action or the running effective 
action. In the above transformation, J{x) is fixed by the 
condition that ^(x) = {ip(x))k and the average is taken 
by using the modified action S'fe[ 7 ’]. 

The running effective action F^)^] continuously inter¬ 
polates between the bare action at the UV scale^^ and the 
exact effective action (or Gibbs free-energy functional) 
F[0], which is the generating functional of the IPI cor¬ 
relation function, when A: —)■ 0. Its evolution with k is 
described by an exact RG flow equation^® 


dTk[c(>] 

dk 


Rk{x-y) (rf 


J xy 


(35) 


with the initial condition rA[/'] = (33[<f)\ and 

T^^\xi, ... ,Xn) = 5'^Tk/5((){xi).. .6(j){xn)- By differ¬ 
entiation, this functional flow equation is equivalent to 
an infinite hierarchy of coupled flow equations for the 
running effective potential, Uk{((i) = Tk[(()]/L^^ and the 
running IPI correlation functions (then all evaluated for 
uniform held conhgurations). 

Finding the exact solution of the functional integro- 
differential equation in Eq. (35) is an impossible task 
in general and one needs to develop approximation, or 
closure, schemes that basically replace Eq. (35) by a 
hnite set of coupled equations for functions. This has 
been systematically and successfully pursued for a series 
of problems in both high- and low-energy physics.® 

Applications of the NPRG formalism to the one¬ 
dimensional theory with S'):/?] given by Eq. (1) have 
been previously considered.In these studies it was 
found that simple approximation schemes fail to recover 
the low-temperature physics of the model, in particular 
the activated scaling of the correlation length. Here, we 
will show what is the underlying reason for this failure. 
To this end, we will hrst derive the exact asymptotic low- 
temperature form of the running effective action F k [/’] by 
using the instanton approach and a mapping to the one¬ 
dimensional Ising model. 


A. Running effective action and instantons in the 
limit T —>■ 0 (C 00) 

We first compute the expression of the running effec¬ 
tive potential, Uk{((>) = ^k{4>\/L with (j) a uniform field. 
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by using the instanton technique in the low-temperature 
regime where the correlation length is large (see also 
above). 

For the one-dimensional theory under study, the 
fc-dependent regularized action reads 


-S'fcM = / dx -{dx<p) +V{‘p{x)) 
Jo 


dx / dy (l){x)Rk{x - y)(j){y) (36) 


- L 


v{M + ^Rkmik 

where V{ip) is given in Eq. (2), 

0 o,fe = argmin^ 


V{v) + 


(37) 


one so that aL = 2n, we get 


poo pioo poo / 

Pk{(t>)=N' dae-“^^^^ / dpdv 

•jQ J— ioo •/0 \i —1 


exp 


y,L{l P aak) + Lv- -W [(/i -b 0^2*-! + (m - 0-^2*] 

(39) 


where M' is another normalization constant. Integration 
over the variables Zi then leads to 


poo pioo 

Pk{(j))=J\f' dae"“^^'^'= / d^dr 

Jo J —ioo 


exp 


L[ fi + afzak + V- -- ln(^^ - 

00 ,fc 2 


(40) 


For large L we can use a saddle point evaluation of the 
integrals over /x and v, which gives 


and we have added the last term in Eq. (36) for conve¬ 
nience, so that S'fe[0o,fe] = 0- Contrary to the previous 
section on the finite-size effective potential, we take here 
the thermodynamic limit and let L —)■ oo. The restriction 
to the spatial extent of the fluctuations is now provided 
by the IR regulator Rk{q)- 

A very simple regulator is the Callan-Symanzik one, 
Rk(q) = which amounts to adding a conventional 
mass term to the bare action. (Note that in this case the 
running effective action is equal to the bare action only 
in the limit A —>■ oo but this has no consequences for 
the physics at intermediate and small momentum scales.) 
It is then easy to see that there exists a threshold k = 
kc such that 0o,fc 7^ 0 for all k < kc- This threshold 
corresponds to the moment along the RG flow where the 
running modified potential V{lp) + develops two 

minima and has a double-well shape. One can expect 
that this qualitative evolution is very general and does 
not depend on the details of the regulator. The precise 
form of Rk{q) changes only the point kc where the double¬ 
well shape first appears. 

For k < kc we can thus evaluate the probability of 
finding a particular magnetization in the system by using 
the instanton method, much like in section IIB. This 
probability is given by 


Pk{4>)=Af' da exp 


- aL/3Sk 


-|-T(a — alna-b — ln[(l + aakY 



(41) 


Finally the integral over a can also be evaluated with 
the saddle point method and the value of a at the saddle 
point is found to be 


a = 



(42) 


where the range of magnetizations is limited to 0 < 0o,fe. 
The low-temperature regime corresponds to j5S1^ large 
and the system is then described by a dilute instanton 
gas. The running effective potential Uk is just —l/(/3L) 
times lnPfc(0), to which one subtracts the contribution 
of the IR regulator, and it is given by 


t^fe(0) = - 4>l,k) + V"(0o,fe) 

= - ^i?fc(0) (0^ - (j)l^k) + ^(00,/c), 

(43) 


Pfc(0 = M/L) = 

n>l 




which is valid for 0 < 0o,fc and is asymptotically exact 
when the temperature goes to zero. The associated flow 
equation reads 


2n 


6 j Zi — L — 2nak j S 


{Z2i-1 — h2i) — — - 


(38) 


where Af is a normalization constant, S'J is the action 
evaluated on an single instanton profile and ak is the in¬ 
stanton width. By exponentiating the Dirac delta func¬ 
tions and passing from a discrete sum to a continuum 


dUkW ^ d_ 
dk dk 

Id 


To go further and find the asymptotic low-temperature 
expressions for the IPI vertices at scale k, we use a short¬ 
cut provided by an approximate mapping when T —)■ 0 





62 




(44) 


K, - 



13 


between the (p'^ theory and the Ising model. The latter 
is described by the Hamiltonian 

L L 

H[{(yi}] = -J^ o-iCTi+i (45) 

i=l i=l 

where ctl+i = cri if periodic boundary conditions are 
used and we have set the lattice spacing to one. The ther¬ 
modynamic limit with L —> oo is considered here. This 
calculation for the one-dimensional Ising model is rather 
standard®^ and the details are given in Appendix B. 

From the comparison between the expressions of the 
pair correlation function obtained in the low-temperature 
limit, where the long-distance properties of both theories 
are described in the continuum, i.e., for the Ising case 

m) ~ (I — . (46) 

with the correlation length 

^(m) ~ i\/l — , (47) 

where m is the magnetization per site, and for the (mod- 
ihed) theory (from the instanton calculation), 

where 


in complete agreement with Eq. (43). 

By using the mapping with the one-dimensional Ising 
model, we can also obtain low-temperature expressions 
for the higher-order IPI correlation functions, 

• • •, at the scale k. With the results given in Appendix C 
we obtain 

rfHpi,P2,P3;0) = (27r)d(pi +7*2 

X [3 - {PiP2 + PiPs + P2P3)] 

^i^\pi,P2,P3,Pi](t)) = (27r)d(pi -fp 2 +P3 

X {-[c{(pf + s{(l)f]^{(j))‘^PlP2P3P4 - [3c{(j)f + s(^)^]x 
i{4'T{PlP2 + P1P3 + PlP4 + P2P3 + P2P4 + P3P4) 
-l- 3 [ 5 c(^)^ -I- s{(j))^]) 

(56) 

where c{(j)) = (j) and c(^)^ -I- s((/))^ = ^0 k should not 
be confused with the notation also used for the pref¬ 
actor of the derivative term in the bare action). More 
generally the running IPI correlation function can be 
cast in the form r^"^(pi, • • • ,p„; </>) = {2tt)6{pi -I- • • • -|- 
Pn)^{4')~^gn{^i4>)pi, - ■ ■ ,i{(l>)Pn',(i>) where the remaining 
(j) dependence in gn does not contain exponential terms 
involving exp(/35^). 

B. Approximation schemes 


one can see that the two theories map onto each other 
with the following formal replacements when T —)■ 0: 

2 J SI (1 — m^) —>■ (<('o,fc ~ 4’“^) (50) 

and the instanton action, is identihed with the 
domain-wall energy when T —)■ 0. 

At low temperature, the two-point IPI correlation 
function can be written as 

^ [p"+r"(3)] - nm 

(51) 

where the last term is due to the definition of the running 
effective action in Eq. (34). This expression can be put 
in the form 

/3-irf (p; <^) = Zk{4)p^ + U'a4) (52) 

where Zi,{(j)) can be obtained as 


Having obtained the exact expressions for the running 
effective potential and the running IPI correlation func¬ 
tions in the low-temperature limit, we can now test ap¬ 
proximation schemes for the exact NPRG equation in Eq. 
(35). Among the several approximation schemes so far 
proposed, we will focus first on the most popular one, the 
so-called derivative expansion. In this approximation, 
the running effective action at the scale k is expanded in 
gradients of the Held, 


^k[4>] = P 


Uki^ix)) + l-Zki^ix)) 


d(j){x) 

dx 


(57) 

where the higher-order terms involve 4, 6, etc., deriva¬ 
tives of the field. We will show that finite truncations 
of the derivative expansion are unable to reproduce the 
exact features of the low-temperature physics. 


Zki4) = lim 4^rf^(p, (l>) = 


1 


p'-T6 2/?9p2-'' 2/3((/)2^fc - (^2) 

In addition, we can check that 


?(</>)■ 




= -/3Sfc 


P 4>lk 4V4lk ’ 


(53) 

(54) 


1. LPA 

The Local potential Approximation (LPA) is the low¬ 
est order of the derivative equation. It corresponds to 


TfcM = /3 


TT t A! \\ ^ ( d4>{x)^ ^ 

Uk{4{x)) + 2 


(58) 
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where the coefficient of the gradient term is constant and tion for values of T higher than the energy barrier of the 

is not renormalized. Plugging this ansatz into Eq. (35), double well, or more precisely than the instanton energy 

computing it for a uniform field (j) and choosing the simple cost S* , they fail to reproduce the low-temperature result 

regulator Rk{q) = bkk^ with bk constant (taken to 1) with a thermally activated dependence of the correlation 

lead to the following differential equation for the running length, oc exp(/35'*). For instance, the curvature of the 

effective potential Uk{4>)'- effective potential in zero, Kfc=o = C^fcLo(O)) which should 

vanish exponentially when T —)■ 0 as exp(—/3S'*) (see 
1 2 n 1 s-lso section II) is generically found to vanish as a power 

dkUk{(l>) — _|_ j ^2 ■ law of T instead. The nonperturbative regime associated 

with the rare localized events, which is captured by the 

.r. 1 , , . . instanton calculation, is therefore completely missed. 

It IS easily verifaed that the exact expression m Eq. (43) 

does not satisfy the above equation. The latter is actually 

unable to reproduce the correct scaling of the correlation ^ derivation expansion 

length, with, e.g., U'f! + k^ cc exp(—/SS”^) [see Eq. (54)]. 

running effective po- corresponds to the following ansatz 

tential Uk{(f>) at several values ol k, as obtained Irom 

the LPA with a regulator of the form"^^ Rk{p) = {k"^ — ^ f ^ 

p^)0(A:^—p^). The curves illustrate the return to convex- P^,[())]=/3 / Uk{(j)ix))-\ —Zfc(^) ( — | . (60) 

ity of the potential. However, as also known from previ- ^ J _ 

ous attempts,if the LPA provides a good descrip- 

When inserted in the exact RG flow equation, this ansatz leads to two coupled differential equations for the functions 
Uk{4>) and Zk{4>) [the latter is obtained from the exact flow equation for the second vertex P), with the use of the 
prescription given in the first equality of Eq. (53)]: 

dkUk{<t^) = Y^dkihkk^) [Zk{rt>) (£/"(</)) + (61) 

dkZk{rt>) = -^dkibkk^) -yj'^'{<j^fZk{cl>f/^ + bkk^y"^ + ^4(().)C/r(0)Zfc(0)-i/2 ([/"(0) + hkk^'^''^ 

+ ([/"(<).) + bkk^Y^'^ - \z'yel>)Zk{^r^l^ {Um + hkk^Y^''^ 

(62) 

where the IR cutoff function is of the same form as above (and a residual /c-dependence is allowed in bk). 

When inserting the exact expression for Uk{(j)) and Zkie)) given in Eqs. (43), (53), and (54), one can see that Eq. 
(61) is now satisfied at leading order in exp(—/35'J) but not Eq. (62). The exact expressions indeed generate terms 
of order exp(2/3S'J) in the right-hand side of Eq. (62) which do not cancel and have no counterparts in the left-hand 
side [which is itself essentially of order exp(/3S'^)]. The problem found at the LPA level can be formally cured at the 
level of the effective average potential but at the expense of an inconsistency at the level of the function Zk{4>). 

3. Fourth order of the derivative expansion 

To check whether the results found above correspond to a more systematic pattern, we have considered the fourth 
order, which corresponds to taking 

Pfc[0] = y Uk{(j){x))+ ^Zk{<j){x)){d(j){xY + ^Yk{(j){x)) {d<p{xY . (63) 

The equation for the running effective potential in Eq. (61) is unchanged but that for Zk{4>) is now obtained as 

dkZk{<)) = ^dkibkk^) ~YJ'Yf>fZk{fYi^ {uY<p) + bkk^Y'^^ + ^z'M)u'Y<l>)Zk{f>Y^i^ {u'YY + bk^Y'^" 
+LzYe)fZk{<i>r^^^ {u'YY + bkYY^'^ - ^z'YYZkYY^i^ {u'YY + hkYY""^ - \YkY)ZkYY^'^^ 
YYY + bkYY^'^ ■ 


( 64 ) 
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An equation for Yk{4i) is also derived by considering the 
flow of the 4-point IPI vertex but it is too long to be 
given here. 

When inserting the exact low-temperature expressions 
for [/fc(</>), (</>), and Yk{(j)) [the latter can be obtained 

from Eqs. (55,56)] in the three flow equations corre¬ 
sponding to the present ansatz, one finds that both the 
equation for Uk and that for Zk in Eq. (64) are sat¬ 
isfied. For the latter, the term involving Yk{(j)) in the 
right-hand side of Eq. (64) now exactly cancels the term 
in exp(2/3S'J) that led to an inconsistency in the second- 
order approximation (see above). On the other hand, one 
can check that the approximate equation for Yk{4>) is not 
satisfied by the exact expression because of the presence 
of terms of order exp(4/3S'^) in the right-hand side [while 
Yfe itself behaves as exp(3/3S'^)]. 


4- General scheme and further approximations 


Guided by the above results, it is now easy to infer 
the general pattern. The prefactors of the terms with 21 
derivatives of the field in the derivative expansion of the 
running effective action rfe[^] are of order exp[(2l—1)/3S'^] 
in the low-temperature regime [and C/^' {(p) + Rk{0) is it¬ 
self of order exp(—,05'^)]. This dominant behavior when 
/3Sl —>■ oo emerges from the exact NPRG hierarchy of 
equations for the IPI vertices because terms that would 
naively lead to a higher power in exp(/3S'^) in the right- 
hand side of the equations (the “beta functions”) exactly 
cancel out. This cancelation effect is however lost if one 
truncates the expansion, whatever the order of the trun¬ 
cation. We conjecture that the appropriate ansatz of 
P k [<p] that reproduces the low-temperature physics of the 
model is instead 


TfcM = 




Yk,2i{(l){x)){d(j){x)) 


21 


( 65 ) 

with, to make contact with the previous notations, 
Yk^2{4>) = and Yk,A{4>) = Yk{(j)). Note that the 

above form of Pfc is not the most general one: in the 
derivative expansion, the term of order is actually a 
combination of terms involving (5(/))^*, 8“^, • • •, 
which even after integration by part cannot in gen¬ 
eral be reduced to a single contribution as in Eq. (65). 
The specific form in Eq. (65) results from the rather sim¬ 
ple momentum dependence of the IPI correlation func¬ 
tions in the one-dimensional Ising model and theory 
at low temperature. 


The above finding allows us to discuss another approx¬ 
imation of the NPRG called BMW.^® It corresponds to 
a closure of the exact NPRG hierarchy at the level of the 


equation for the IPI two-point function P), ’{p,(p): 

<P) = PJ ^dkRk{q) [Gkiq,(p)‘^Gk{p + q, (p)x 

rf ^ (p, q, -p - q; 4>f - \Gk{q, (p, -p, q, -q; (p)] 

( 66 ) 

where Gk{p,(p) = [T^^\p,(p) + Rk{p)]~^ ■ The BMW 
closure consists in setting to zero the internal mo¬ 
mentum q appearing in the 3- and 4- point ver¬ 
tices in the right-hand side. After using the consis¬ 
tency relations, r®(p, 0, —p;^) = dT^f^\p-, cp) / dcp and 
r‘'j^\p,—p,0,0;(p) = d‘^T‘'^\p; (p) / dep"^ , one obtains a 
closed equation 

dk^^k\p^ = ^dkRk{q) [Gk{q,(p)‘^Gk{p + q, (p)x 

dr^^\p;(P)'\ 1 ,.2d‘^^k\p;(p)i 

d<P \ -2^'=^'?’'^) d<P^ J 

(67) 

which can be combined with the equation for the running 
effective potential Uk{(p)- It is easily checked that Eq. 
(67) is not compatible with the exact low-temperature 
expressions of Uk and P^^^ given in section III A: after 
scaling the momenta by ^ {(p) (see section III A), the left- 
hand side of Eq. (67) scales as whereas the right- 
hand side has a term in that does not cancel out. Just 
like truncations of the derivative expansion, the BMW 
closure is therefore unable to properly describe the non- 
perturbative physics of the one-dimensional at low 
temperature. 

The alternative to the existing approximation schemes 
of the NPRG is to start from the exact low-temperature 
ansatz in Eq. (65). This however leads to an infinite 
set of differential equations that cannot be treated with 
standard methods. We have tried another route which 
amounts to considering the running effective action as 
being local in the two variables (p{x) and d(p{x) and intro¬ 
duce an auxiliary field (p{x) to decouple d<p{x) from <p(x). 
This procedure however is highly ambiguous. In addition, 
say we end up with a running effective action of the form 
Tk[(p,(p] = J^Vk{(p{x),(p{x)), it is not clear that stan¬ 
dard approximations on this ansatz will correctly capture 
the expected low-temperature physics. Actually we have 
tried an LPA approximation at the level of the two fields 
(p and (p and it completely misses the nonperturbative 
regime. More work is needed to possibly find a solution 
to this unsatisfactory theoretical situation. 


IV. DISCUSSION AND CONCLUDING 
REMARKS 

In this work we have studied the theory at low tem¬ 
perature in the regime where the behavior of the system 
is completely dominated by nonperturbative instantonic 



16 


fluctuations. We have first discussed empirical recipes 
to extract the parameters of the underlying microscopic 
or effective theory from the numerical study of finite-size 
systems. This strategy could be very useful in the anal¬ 
ysis of finite-size numerical simulation of glassy systems 
but the application to this problem deserves further work. 

Our study also illustrates the difficulty to describe 
the low-temperature nonperturbative physics of the one¬ 
dimensional theory through truncations of the NPRG. 
In a sense, however, the one-dimensional case is harder 
than the situation in higher dimensions. There, the tran¬ 
sition associated with a spontaneous symmetry breaking 
is not destroyed by the fluctuations and the return to con¬ 
vexity has been shown to be properly described through 
simple approximations of the We now dis¬ 

cuss in more detail this higher-dimensional situation. 

Consider for instance the 3-dimensional theory. As 
far as the finite-size effective potential Ul{(I>) is con¬ 
cerned, one can repeat and adapt the qualitative argu¬ 
ments developed in section II A. At low enough temper¬ 
ature, the bare potential has two minima in, say, (j) = ±I 
and the relevant excitations above the uniform ground 
states are system-spanning domain walls or interfaces be¬ 
tween regions of essentially constant positive and nega¬ 
tive magnetization. When the system size L becomes 
larger than the interface width, the system can accommo¬ 
date one system-spanning interface: Ul{4>) should then 
have, on top of the two symmetric minima for ~ ± 1 , 
a plateau for intermediate values of the field; the height 
of the plateau compared to the bottom of the minima 
is given by T /L where T is the surface tension. As L 
increases, this height decreases and goes to zero in the 
thermodynamic limit. The effective potential is convex 
with a flat intermediate portion corresponding to phase 
coexistence. The evolution with L of UL{(j>) is schemati¬ 
cally depicted in Fig. II. In this case, studying finite-size 
systems should allow one to extract two physical quanti¬ 
ties, the surface tension and the correlation length which 
corresponds to the interface width. 

This 3-dimensional theory in the symmetry-broken 
region®® has also been studied in detail within the NPRG 
framework.Although influenced by domain walls as 
in one dimension, the long-distance physics is nonethe¬ 
less different as these nonperturbative fluctuations are 
not strong enough to destroy the phase transition. As a 
result, simple approximation schemes of the running ef¬ 
fective action properly capture the effect of these fluctua¬ 
tions. The running effective potential Uk{(t>) evolves with 
decreasing k from the bare double-well potential to a con¬ 
vex effective potential when fc = 0 : see the schematic plot 
in Fig. II. Provided one chooses an appropriate class of 
IR regulator,® the intermediate “inner” part of Uk{<p) dis¬ 
plays at small k a parabolic shape oc k'^cj)'^ that comes in 
addition to the two symmetric minima in ^ = ±^o,fe- 
This parabolic dependence corresponds to the expected 
exact behavior obtained by considering the nonuniform 
configurations of the field involving domain walls. The 
remarkable feature is that this behavior is recovered by 



FIG. 11. Schematic plot of the evolution of the shape of the 
finite-size effective potential 17 l( 0) (top) and of the running 
effective potential Ukitj)) (bottom) with either system size L 
or running IR momentum scale k for the theory in 3 di¬ 
mensions. 

using approximations of the NPRG, such as the first or¬ 
ders of the derivative expansion, which only consider ex¬ 
pansions about uniform fields.^® This is in stark contrast 
with the situation encountered in one dimension. 

The nature of the NPRG flow somehow changes when 
the running IR momentum scale k crosses some value k^, 
that roughly corresponds to the point at which k^ be¬ 
comes of the order of magnitude of the curvature of the 
running effective potential in (/> = 0: kl (0)|; this in 

turn corresponds to the point where I/fc becomes of the 
order of the width of the domain wall in the (nonuniform) 
field configurations that minimize the running effective 
action at this scale Whereas such information can 

be included in an improved instantonic theory of nucle- 
ation in the cases where metastability in present,®’®® as, 
e.g., when applying a nonzero external source or mag¬ 
netic field, its interpretation in terms of physical quan¬ 
tities of the actual, macroscopic system remains unclear. 
In particular this length scale 1/k* does not provide any 
direct information on one of the important length scales 
in nucleation problems, be., the size of the critical droplet 
or bubble. 

In any case, we think that providing a generic solution 
to the problem posed by nonperturbative fluctuations in 
model systems as the one studied here would be very 
profitable for tackling the harder situations encountered 
in glassy systems which involve activated dynamics in a 
complex landscape. 
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Appendix A: Computation of the combinatorial 
factors for the gas of instantons 

The combinatorial coefficients l 2 n{L) are configuration 
integrals of 2n domain walls of width cr on a ring of size 
L. This problem is equivalent to the computation of the 


partition function of 2n (discernible) hard spheres of size 
CT on a ring of size L m. D = 1. As done in the main 
text, we define Xi, i = l,...,2 n, as the lengths of the 
regions with constant <j) = ±1 {i.e., the gaps between 
the spheres). These variables must satisfy the constraint 
+ 2ncr = L. 


One then has 


hn{L) 


pL—2n(7 

rL—2n<7—X2n-i 

l‘L—2n<7 — (x2n-l+X-in--iA - \-X2) 

/ dX2ri-l 

/ da;2„-2 • • 

• / dxi 

'0 

'0 

Jo 


where the factor L comes from translational invariance and the factor 1/n accounts for the number of ways one can 
choose the first kink/anti-kink pair. In the following, we will determine the expression of l 2 n by recurrence. In order 
to do this, it is convenient to introduce the functions 


9n{y) 


ry 

1 

3 

1 

|•y-iXrl-l+X„-2-\ - 1-2:2) 

/ dXn-l 

/ da;„_2 • • 

• / dxi 


/o 

Jo 


(Al) 


in terms of which the combinatorial factors can be ex¬ 
pressed as 

hniL) = - g 2 n{L - 2na). (A2) 

n 


As a result, the integrals over the variables Xi can be di¬ 
vided into separate integrations over even and odd gaps, 
which can be written in terms of the functions gn(jj) de¬ 
fined above. This yields 


From the definition in Eq. (Al), one can write + 1 ( 2 /) 
in terms of g„(y), 

ry ry 

gn+l{x) = / dXngn{y-X„)= / dXn gn{x„) , (A3) 

Jo Jo 

which, from ( 72 ( 2 /) = Jq dxi = y and by recurrence, im¬ 
mediately leads to 


J2n{M, ~ ^9n 


L — 2na — M 
2 

L — 2na + M 
2 


X 9 . 


»"'!') = (IvTjT ■ (A'l) 

Finally, after plugging Eq. (A4) into (A2), one obtains 
Eq. (13) of the main text. 

In order to compute the combinatorial factors 
<^ 2 n(Af, L) one has to impose that the gaps Xi satisfy 
the two following constraints: 


The extra 1/2 factor comes from the fact that only half 
of the configurations, namely, those with the first domain 
wall joining ip = —1 to ip = -1-1, contribute to magneti¬ 
zation +M, whereas the others contribute to —M. After 
using the exact expression in Eq. (Al), one finally finds 


' 2n 

Xi = L — 2na , 

i=l 

n 

(a;2*-i - X2i) = M, 

< 2=1 


J2„(M, L) 


L 

2n 


{L - 2naf - 
22 ("-i)(n - 1)!2 


(A5) 


which can be rewritten as 


'^X2i-l = 


L — 2ka + M 


y]a:2^ = 


L-2ka- M 


which, with the help of the intensive variables (j) = M/L 
and a = ajL, leads to Eq. (18) of the main text. 
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In the following we show that PLi4>) given in Eq. (18) is properly normalized to 1. We start by computing the 
integrals over (p of the terms of the sum separately. By changing variable to x = (p/{I — 2na) one gets 


T2u{L) = 2 


2n ^ + (1 —2na) 


-(l-2na) 




(1 — 2 na)^ — (p"^ 


n!(n — 1 )! 


= 2 


(C/2) 


2n 


n!(n — 1 ) 


^ + 1 

, (l-2na)2("-i)(l-2na) J dx)!-^^) 


2\n-l 


The integral over dx can be computed as 

r + l 


^ + 1 

J dx(l-x2)"-i =2y dd(cosd)2”-i = Vtt 


r(n) 


r(n + l/ 2 ) 


where r(n + 1/2) = 2 "(2n — 1)!! -y/Tr. By using the fact that 

( 2 n)! ( 2 n)! 


( 2 n- 1 )!! = 


( 2 n)!! 2 "n! ’ 


one then hnds 


T 2 „(T) = 2 C 

From Eqs. (18), (14) and (A 6 ), one ends up with 

1 


2 „ (1 - 2na) 


2n-l 


( 2 n)! 


r+l 


/ dpPLiP)= . . 

/-I ZL{C,a) 


1 


„ + l l/( 2 a) 

J (i(() (d((() — 1) + (5((/) + 1)) + ^ 2^ 


-1 

1/C2a) 


ZL{C,a) 


2(^2n (1 - 2na)^" ^ ^ ^ 


n =0 


( 2 r 


2 „ (1 - 2na) 


2n-l 


( 2 n)! 


(A 6 ) 


(A7) 


The expression of Pl{(P) for th® one-dimensional Ising 
model^^ can be recovered as a particular case of Eq. (18) 
in the limit a 0 {i.e., for infinitely sharp domain walls) 
and for S* = 2 J. In particular one then has 


ZL{C,a = 0) = 2j2 


e 


n—O 


( 2 n)! 


= 2 cosh ( , 


and 


PLi<l>) = 


1 


2 cosh ^ 


(j(^_l) + ^(^+l)) 

(l-(/)2)"-l 
n!(n — 1 )! 


+ E 2 (C/ 2 )^ 


These expressions coincide with the results of Ref. [45]. 


Appendix B: Instanton calculation for the 
one-dimensional Ising model 

As stated in the main text, we consider the one¬ 
dimensional Ising model with periodic boundary condi¬ 
tion which is described by the Hamiltonian 

L L 

H[{ai}] =-j'^aiai+i - h'^ai (BI) 


where ctl+i = ai, the lattice spacing which is as unity, 
and, contrary to the case studied in section II, only the 
thermodynamic limit L —>■ oo is considered. 

We summarize the main (known) results about the 
model. The partition function can be computed using the 
transfer matrix method.®"^ The transfer matrix is given 

by 


V = 


'^h+j g-j ' 

g-J (,-h+J 


(B2) 


where h = (3h and J = f3J and the partition function is 
obtained as Zm = TrY-^. One can easily diagonalize the 
transfer matrix with the rotation 


U = 


/ cos 9 — sin d 
I sin 9 cos 9 


(B3) 


where 1/tan 2d = e^'^sinh/i. The eigenvalues are given 
by 


A± = cosh h ± Y cosh^ h — 2 sinh 2 J . (B4) 

The average magnetization m =< Oi > is then 

sinh h 


m = 


\J cosh^ h — 2 sinh 2J 


(B5) 
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and the two-point connected correlation function where 

= sm^{2e) , (B6) 

which can be rewritten as 



(C3) 


=sm‘^{29)e-G^ (B7) 

-1 

. These expressions provide the 
magnetization and the two-point function at fixed exter¬ 
nal magnetic field h. However we would like to have the 
magnetization instead of the external field as the primary 
variable since we want to work with the effective action 
defined by the Legendre transform. (34) 

One can thus invert the relation in Eq. (B5) to obtain 


where ^ = In f 


After diagonalizing the transfer matrix (see Appendix B) 
and using that 


U-I 5 U 


/ cos 26* — sin26*\ 

sin 29 — cos 29J ’ 


(C4) 


one obtains in the thermodynamic limit 


A+ 

A_ 


W 1 — m‘^(l — + e 


(B8) 


1 — m^(l — — e 


-2J 


so that, in the limit of very low temperature, one gets 

2 e- 2 ^ 


^(m) ~ 


VT^ 




(B9) 


and 


G(2) (r) = (1 - . (BIO) 

These results are used in section III A. 


Appendix C: Derivation of and F))*^ from the 
one-dimensional Ising model 


d4) 




ri+r2'' 


(C5) 


where c = cos 29 = m and s = sin 29 = 1 — (c should 
not be confused with the notation also used for the pref¬ 
actor of the derivative term in the bare action). The 
connected 3-point correlation function is then given by 


X (<7i-t-ri_+r2 ~ (®’i+i’i+r 2 ))) 

^ (t) ■ 

(C6) 


We start from the calculation of the 3-point correlation 
function in the Ising model. We need to compute 

(o'iO'i-j-riO'z-}-r^+r2) • (^4) 

It is given by 

(a.a.+.,a.+.,+.,) = Tiv^Tr 

A_|_ -T A_ 

(C2) 


If we call xi = i, X 2 = i + ri, X 3 = i + ri + r 2 , the above 
result translates into 

/A \4^3“4:i 

W^^^(xi<X 2 <X 3 ) = -2cs^(j^] , (C7) 

so that the 3-point connected correlation function for 
generic arguments can be written as 


W^{xi,X2jX3) = - 2 cs^ 


0(a;3 - a;2)0(a;2 - xi)-I- f 0(a;i - a;2)0(a;2 - xa) 


"A /\ 

+ ( ^ j 0(a^2 - a:3)0(a:3 - a;i)-b f ^ 


0 (a;i - X3)0(a;3 - X 2 ) 


(C8) 


' \ ai2— / y^ \0C3—X2 

+ ( ^j 0(a;2 - a;i)0(a;i - a^a) + f 0(a;3 - a;i)0(a;i - 0 : 2 ) 
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where 0(a;) is the Heaviside step function. Neglecting the underlying lattice and performing the Fourier transform 
lead to 


W^^'>{pi,P2,P3) 


(27r)(5(pi +P 2 +P 3 ) 


4cs 2^ ^ ((piP 2 +P 1 P 3 +P 2 P 3 ) - 3C 


(C9) 


We now use the mapping between the Ising model and the (p^ theory at low temperature and the relation between 
the connected 3-point correlation function and the IPI 3-point vertex.^® We finally obtain 


rf^(pi,P2,P3) = -ff (Pl)f(2^(p2)ff (P3 )W^^'Hp1,P2,P3) 

= (27r)(5(pi +P2+ Ps) 

X - (piP 2 +P 1 P 3 +P 2 P 3 )] 


where c = (j> and s® = — (jP. Note that it is 

obtained from = rfc[(/)] -b 

I f dqRk(q)(l>(g)(l>(—q), which appears in Eq. (CIO) and 
not rj,^\p). The calculation of F^"*^ can be done in an 


analogous way and leads to Eq. (56). Although a cum¬ 
bersome derivation, the higher orders can also be ob¬ 
tained along the same lines. 
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